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ABSTRACT 

Context. Hard X-rays from solar flares are an important diagnostic of particle acceleration and transport in the solar atmosphere. 
However, any observed X-ray flux from on-disc sources is composed of direct emission plus Compton backscattered photons (albedo). 
This affects both the observed spectra and images as well as the physical quantities derived from them such as the spatial and spectral 
distributions of accelerated electrons or characteristics of the solar atmosphere (e.g. density). 

Aims. We propose a new indirect method to measure albedo and to infer the directivity of X-rays in imaging using RHESSI data. We 
describe this method and demonstrate its application to a compact disc event observed with RHESSI. 

Methods. Visibility forward fitting is used to determine the size (second moment) of a disc event observed by RHESSI as a function 
of energy. Using a Monte Carlo simulation code of photon transport in the chromosphere, maps for different degrees of downward 
directivity and true source sizes are computed. The resulting sizes from the simulated maps are compared with the sizes from the 
observations to find limits on the true source size and the directivity. 

Results. The observed full width half maximum (FWHM) of the source varies in size between 7.4 arcsec and 9.1 arcsec with the 
maximum between 30 and 40 keV. Such behaviour is expected in the presence of albedo and is found in the simulations. The uncer- 
tainties in the data are not small enough to make unambiguous statements about the true source size and the directivity simultaneously. 
However, a source size smaller than 6 arcsec is improbable for modest directivities and the true source size is likely to be around 7 
arcsec for small directivities. 

Conclusions. While it is difficult to image the albedo patch directly, the effect of backscattered photons on the observed source size 
can be estimated. This is demonstrated here on observations for the first time. The increase in source size caused by albedo has to be 
accounted for when computing physical quantities that include the size as a parameter such as flare energetics. At the same time, the 
study of the albedo signature provides vital information about the directivity of X-rays and related electrons. 
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1. Introduction 

X-ray emission from solar flares is an important diagnostic 
of electron acceleration and transport in the solar atmosphere. 
In the common flare scenario, particles are accelerated in the 
corona. They precipitate along the field lines of a magnetic loop 
to the dense chromosphere where they loose their energy in 
Coulomb collisions with the surrounding plasma and produce 
hard X-ray (HXR) bremsstrahlung radiation. However, the de- 
tailed physics of electron acceleration and propagation is a sub- 
ject of active research. Many studies of the different aspects of 
the physics involve the determination of the spatial characteris- 
tics (position, size) over which processes in flares happen. An 
accurate measurement of the size of X-ray flare sources is cru- 
cial for the computatio n of quantities such as the total energy 
involved in a flare ( eg. Canfield et al. 1980; Strong et al. 1984; 
Saint-Hilaire & Benz(l2005t IVeronig et alj|2005bt lHannah et alj 
2008; Fivian et al. 2009 ) or densities of coronal sources (eg. 
Batta glia et alJl2009l:lKrucker & Linll2008l) . 



The Reuven Ramaty High E nergy Solar Spectroscopic 
Imager (RHESSI. iLin et aTll2002l) is capable of observing X- 
rays from solar flares with high spectral and spatial resolution. 
This allows us to determine the positions and sizes of HXR flare 
sources to very high accuracy. lAschwa nden et al. (2002) mea- 
sured the position of X-ray sources as a function of energy, find- 
ing a decrease of the radial position and size with increasing en- 
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ergy. iKontar et"aD (120081) used the new analysis method of vis- 
ibility forward fitting dHurford et alJl200%ISchmahl e t al. 2007) 
on a limb event, finding a clear decrease of radial position with 
energy. It is now possible to not only find the position but also the 
sizes of the sources par allel and perpendicu lar to the magnetic 
field as demonstrated bv lKontar et al.l ( 120101) with visibility for- 
ward fitting. The parallel sizes provide valuable information on 
the physics of electron transport in the chromosphere. In a target 
with steeply increasing density in which the electrons predomi- 
nantly lose their energy in Coulomb collisions, the parallel size 
is expected to decrease with energy. The size perpendicular to 
the field lines contains information about the convergence of the 
magnetic field. Since higher energetic electrons penetrate deeper 
into the chromosphere, a decrease of the perpendicular size with 
energy is also expected for a converging m agnetic field. Such 
behaviour was found by Kontar etalJ i fcoiOh in a RHESSI limb 
flare. 

RHE SSFs highest an gular resolution in theory is 2.3 arc- 
sec dHurford et alJ l2002h . However, typically observed source 
sizes (in terms of FWHM ) are often in the ra nge of 10 arc- 
sec fo r coronal sources (eg . lKrucker & Linl20 08) and microflare 
loops (jHannah et al]l2008l) and around 3-7 arcsec for footpoints 
dDenn is & Pernak 2009). In comparison, measurements of flare 
ribbons in EUV (eg.lWang et al.ll2007l:lFletcher & Warrenll2003h 
and white light ( Iso be et al.l2007b " Hudson et al.l2006l) find much 
smaller structures on the scales of 1 arcsec. One common ex- 
planation for this is that RHESSI might not fully resolve those 
smaller structures. 
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A more likely explanation is X-ray albedo: downward di- 
rected photons can be Compton backscattered in the photo- 
sphere and add to the observed X-ray flux. The influence of 
this on HXR spectra was f i rst investigated by TombU^ d 19721) 
and lSantangelo et al.l dT973b.lBai & Ramatvld 19781) and more re- 
cently IZhang & Huang! (|2004|) and IKontar et all d2006l) showed 
that this albedo flux can account for up to 40 % of the de- 
tected flux in the energy range between 30 and 50 keV even in 
isotropic sources. The backscattered component originates from 
a large area compared to the size scale of typical flare sources. 
Thus the brightness of t his albedo-patch is rather fa int and di- 
rect imaging is difficult (ISchmahl & Hurfordl 1200 2). However, 
IKontar & JeffrevI (120101) recently found how the albedo patch 
will affect the measured size of the source even if direct imag- 
ing is not possible. Using Monte Carlo simulations of photon 
transport they calculated the expected source sizes and positions 
of X-ray sources on the solar disc, taking into account albedo. 
They found that the observed size depends on the location on 
the solar disc, the spectral hardness and the downward directiv- 
ity (ratio of downward to upward photon flux a = Idownlhip)- 
Furthermore, the size as a function of energy increases with in- 
creasing energy, peaking at around 40 keV, before decreasing. 
However, albedo is not just an unwanted effect that has to be cor- 
rected for but can be used to study the directivity of X-ray emis- 
sion in solar flares. From the knowledge of the X-ray directiv- 
ity, information on the directionality of the ac celerated electron 
distrib ution can be gained as demonstrated by dKontar & Brownl 
2006). Currently, this is the only available diagnostic of down- 
ward electron beaming, except for contradictive measurements 
of po larization (see eg. ISuarez-Garcia et aT. 2006; Bogg s et al.l 
2006). Numerical modelling (eg. lLeach & Petrosianll 19831) sug- 
gests values of a in the range of 1-2. Rec ent observations of the 
X-ray directivi ty indicate values of 0.3 -3 (Kaspar ova et al.l2007l) 
and around 1 (Kontar & Brown 2006|). 



None of the previous studies on directivity and electron 
beaming used imaging observations to measure the directivity. 
Since the albedo patch is faint it is difficult to image directly. 
Furthermore it is hard to interpret different source contributions 
correctly in single images and to distinguish potential extended 
primary sources from albedo. In this work we discuss an indi- 
rect method, based on imaging spectroscopy of the moments of 
the HXR spatial distribution. Visibility forward fitting is used 
to determine the source size (second moment of the spatial X- 
ray flux distribution) as a function of energy. Since albedo is 
energy dependent, this should be reflected in the observed size. 
The method is applied to observations of an on-disc flare, as well 
as simulated pho t on ma ps. The Monte Carlo code developed by 
IKontar & JeffrevI (|2010) was used to simulate photon maps for 
different true source sizes and degrees of directivity. From those 
maps, X-ray visibilities were created and fitted using visibility 
forward fitting the same way as the observed data was analysed. 
Since albedo is energy dependent with a peak around 40 keV the 
observed size of the source as a function of energy will follow 
a distinct pattern. Comparison of this observed pattern with pre- 
dictions from the simulations makes it possible to find a measure 
for the true source size and the directivity. 

In Section [2] we describe the new method. Section 12.11 de- 
scribes how simulations of photon transport were performed for 
comparison with the observations. In Section [3] we apply the 
proposed method to observations of an actual RHESSI flare fol- 
lowed by a discussion (Section|4]i. 



2. Method 

For an X-ray source of brightness I{x,y) [photons cm~ 2 
arcsec~ 2 ], the source size can be obtained from the second mo- 
ment of the brightness distribution: 
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where we integrate over the whole emission area S. The source 
size (FWHM) is the n given as FWHMr V _ 2 y/2ln2(x size /y size ). 
In their simulations. IKontar & JeffrevI (12010) computed the mo- 
ments directly from the simulated photon maps. For imaging ob- 
servations of true flares with RHESSI, this is impractical for 
various reasons (background, CLEAN side lobes, artefacts of 
the imaging algorith m). We therefore use the method of visi- 
bility forward fitting dHurford et al.ll2002t ISchmahl et aI1l2007l) 
to determine the FWHM of the source. Previous studies demon- 
strated the usefulness of visibilities to mea sure the moments o f 
the HXR-distribut i on (ze roth moment: flux, Ems lie et al.ld2003l): 
Batta glia & Benzl y()()6). first moment: posi t ion, [K ontar et al 
(20081), second moment: size, IKontar et al.l d2O10l) ; Ixu et al 



(2008)). Visibility forward fitting is also ideally suited to look for 
albedo because it emphasises large scale structures in a source. 
This can be shown as follows. The second moment of the spa- 
tial X-ray distribution as defined in Eq. Q] relates directly to the 
size found in visibility forward fitting. The X-ray visibilities of 
the source are the sp atial Fourier transforms of the photon flux 
dHurford et alJ l2002): 
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or in vector form 
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where k - (u,v),r = (x,y). The total, spatially integrated X- 
ray flux is simply represented by V(k) at k — 0. Indeed V(k = 
0) = L I(r)dr as can be seen from Eqs. [2] and [3] The size of the 
source is directly related to the second order derivative of Eq. [2] 
at k = 0: 
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Using Eqs. |4] and |5] in Eq. Q]the size of me HXR sources from 
X-ray visibilities is: 
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Let us assume a Gaussian source I(x,y) - J±- e ( -~ ( ' x ~ +r} l 2cr ~\ 
where the centroid is at x = 0, y = 0. The size of the source is 
represented by the Gaussian sigma. The Fourier transform of a 
Gaussian is also a Gaussian: V(u, v) ~ e - 27r '°~( ir+v '\ Substituting 
this in Eqs. [6] and [7] one gets the sizes of the HXR sources: 
x 2 - = cr 2 and y 2 . = cr 2 . Thus, fitting Gaussian sources to 
RHESSI visibilities naturally emphasises the large-scale struc- 
tures of any flare source and the resulting size is directly related 
to the second moment of the HXR emission. 
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Fig. 1. Simulated photon map I(x,y) for true source of size 8 arc- 
sec FWHM (white dashed circle) at energies 30-45 keV. The 
white line indicates the FWHM of a circular Gaussian source 
that was forward fitted to visibilities corresponding to the simu- 
lated map. 

2.1. Comparison of observations and modelling 

To find a measure of the true source size and constrain the 
directivity, measurements of the observed FWHM have to 
be compared with pre dictions from the code developed by 
iKontar & Jeffrevl (HblO). The following steps are taken: (i) for a 
given true source size and downward directivity, a photon map of 
the total (direct and backscattered flux) is produced for a source 
at a heliocentric angle corresponding to the location of the ob- 
served flare (see Fig. [TJ. (ii) The photon map is Fourier trans- 
formed to get the corresponding visibilities, (iii) Those visibili- 
ties are fitted using visibility forward fitting to find the FWHM, 
in the same way as for the real data. This size is directly com- 
parable to the FWHM found for the data. The simulated size 
depends on two free parameters, the true source size and the di- 
rectivity. Thus one cannot infer either of those parameters from 
direct fitting. One can find one assuming the other and compare 
the resulting size behaviour to the observations, finding upper 
and lower limits of the directivity and the true source size. 

3. Observations 

In this section we present observations of a compact flare to 
which we applied the above described method. The flare pre- 
sented here was a compact GOES C2.2 class event on the so- 
lar disc. It was observed on the 2nd September 2002 with peak 
time 00:08 UT at a position of E31N6. This corresponds to a 
heliocentric angle of 6 = 30° (or expressed in terms of the co- 
sine: fi = cos(0) = 0.86). The event is well suited for analysis 
by the moment based approach since it has a simple structure. 
Furthermore, the albedo effect is more pronounced close to the 
centre of the solar disc. Thus, a significant albedo contribution 
is expected. Figure [2] (inset) shows the RHESSI lightcurve of 
the event at 6-12 keV and 25-50 keV. Full Sun spectroscopy 




Energy [keV] 



Fig. 2. Spectrum of the observed event, fitted with a thermal 
component (blue) plus a thick target power-law model (red). 
Purple gives the background level. Inset right-hand corner: 
Lightcurve of the event at 6-12 keV (thin line) and 25-50 keV 
(thick line). The dashed lines indicate the observed time inter- 
val. 



was used to determine the shape of the spectrum (needed for 
modelling), as well as the energy range of the thermal and non- 
thermal emission. Figure [2] displays the spectrum fitted with a 
thermal plus thick target model with albedo correction for an 
isotropic source. We find that the spectrum is non-thermal above 
about 16 keV with an electron spectral index <5 = 4.1 correspond- 
ing to a photon spectral index y — 3.1 in a thick target. 

3.1. Imaging 

3.1 .1 . Visibility forward fitting 

Four energy bands (16-20, 20-30, 30-45, 45-75 keV) of the non- 
thermal emission (Fig. |2]i were used for the analysis. The en- 
ergy bands were chosen in order to cover the energies with non- 
thermal emission, yet allowing for a high enough signal to noise 
ratio in each energy band. A time interval of 44 seconds was 
used, which spans the full HXR peak emission (dashed lines in 
the inset in Fig. [2). The source was fitted with a single circular 
Gaussian model. The visibility fits are shown in Fig.[3]for each of 
the four energy bands. The top half of each panel shows the mea- 
sured visibility amplitudes with errors and the corresponding fit 
while the bottom half gives the normalised residuals. The visi- 
bilities suggest negligible modulation in the finest grid (1), and 
only small modulation in the second finest grid (2), suggesting 
a lack of detectable structure at the spatial scales of detectors 1 
and 2. This gives a first lower limit of the source size between 3.9 
and 6.8 arcsec. Fits using all grids compared to using only grids 
3-9 gave equal results within uncertainties, further confirming 
that the contribution of the finer grids is small. Also, the chosen 
number of roll bins (coverage in the uv-plane) does not affect 
the fitted size. While a large number of roll bins might be impor- 
tant for the measurement of higher moments (e.g. shape of the 
source), the uncertainties of the smaller moments (position,size), 
are larger than any effect a different set of roll bins has on the fit- 
ted parameters. 
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Fig. 3. Visibility fits for the four energy bands. The top half in each panel shows the measured visibility amplitudes (black dots) 
with statistical errors (blue lines) and the model fitted to the visibilities (red). The bottom half in each panel gives the normalised 
residuals. 



3.1.2. Other imaging algorithms 

In ad dition to visibi lity forward fitting , CLEAN dHurford et al.l 
I2002I) and Pixon dMetcalf et alj fl996) images were also gener- 
ated for a better understanding of the source structure. Figure [4] 
displays contour plots of the flare at different energy bands us- 
ing these 3 imaging algorithms. Pixon indicates a substructure 
in the 30-45 keV energy band. The extent of the separation of 
the two sources strongly depends on the Pixon-sensitivity (how 
sensitive Pixon reacts to small scale structures) and grids used 
(for the images in Fig. [4] grids 2-8 were used). It could be ar- 
gued that these are two footpoints of a loop but it is not possi- 
ble to make any unambiguous statements. The existence of two 
sources is confirmed to some extent by CLEAN. However, this is 
quite sensitive to the parameters used in the image reproduction. 
The substructure is apparent in CLEAN when uniform weight- 
ing (more emphasis is given to the finer grids) is used instead of 
natural weighting. On the other hand, it appears in images when 
using grids 2-8 with natural weighting and by changing the width 
of the CLEAN beam (a clean beam width factor of 2 was used 
for the generation of the images in Fig . |4j . In the highest energy 
band (45-75 keV), CLEAN is affected by noise and the 30 % 
contour in Fig. [4] outlines some of those artefacts. See Section 
I4. 1 I for a more thorough discussion of how this might affect our 
results. 

3.2. Results 

Figure |5]displays the observed source size found from visibility 
forward fitting (data points). The FWHM was found to increase 
from 8.5±0.4 arcsec in the 16-20 keV energy band to 9.1 + 0.3 
arcsec in the 30-45 keV energy band, then decreases to 7.4 ± 
0.9 arcsec. This was compared to simulations as described in 
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Fig. 5. Observed FWHM as a function of energy (data points 
with error bars). FWHM from visibility forward fitting to sim- 
ulated data for 3 primary source sizes and different degrees of 
directivity (see legend in the figure). 



Section f2TI For the simulations, a heliocentric angle /v between 
0.8 and 0.9 was used and a source height of 1 Mm above the 
photosphere. Due to the large uncertainties in the data, no un- 
ambiguous statements about either directivity or true source size 
can be made. However, one can give a lower limit of the true 
source size, as well as an estimate of the directivity. Two source 
sizes were simulated (6 arcsec and 7 arcsec) for various degrees 
of directivity. The orange, blue and purple curves in Fig. [5] in- 
dicate the expected size for a true source size of 7 arcsec and 
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Fig. 4. Contour maps (30, 50, 70, and 90 % of maximum emission) from three different algorithms (top to bottom: CLEAN, Pixon, 
visibility forward fitting) and energy bands (left to right: 16-20 keV, 20-30 keV, 30-45 keV, 45-75 keV). 



directivities a — 1 , a = 2 and a — 3 respectively. The red curve 
gives the expected size for a true source size of 6 arcsec and di- 
rectivity a — 5. This excludes a very small source since very 
large directivities would be needed to explain the observed size. 
The dashed green curve indicates a model for which the directiv- 
ity is increasing with increasing energy as a — e/10 where e is 
the photon energy in keV. Such a behaviour would be expected in 
a collisional target since electrons with smaller energies are ex- 
pected to have a higher number of collis ions which w ill result in 
a larger spread of the pitch angle (see eg. lBrownll972h . Finally, a 
hypothetical constant source size is indicated by the yellow line. 
Table Q] provides the probabilities that a cho sen model is con - 
sistent with the data according to the^- 2 test dPress et alj|1992l) . 
On these grounds, a true source size of 7 arcsec in the case of 
isotropic albedo can be rejected confidently, as can the case of 
increasing directivity with energy. Previou s observational stud- 
ies found directivities a round a = 1 to 2 (Kasna rova et al.l l2007: 
iKontar & Brownl 12006). Provided that this is a reasonable as- 
sumption for the flare presented here, the true source size is then 
most likely of the order of 7 arcsec for a directivity of a — 2. 



Table 1. Probabilities that a given model is consistent with the 
data according to the^- 2 test. 

Model Probability that model is 

consistent with data 

FWHM=7" a = 1 02% 

FWHM=7" a = 2 56 % 

FWHM=7" a = 3 24 % 

FWHM=6" a = 5 24 % 

FWHM=7"a = e/10 6% 

FWHM=8.8" no albedo 42 % 



4. Discussion 

The method described here used the energy dependence of 
albedo and the ability of visibility forward fitting to get an in- 
direct measure of albedo in imaging. This is an advance com- 
pared to previous attempts in which albedo was only tried to 
be imaged in one energy bands. The limitation of the presented 
method is clearly the large uncertainties in the measured size 
that do not allow an unambiguous fit of any model to the obser- 
vations. However, they do provide important constraints on the 
source size and directivity. Further, the morphology of the anal- 
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ysed source has to be investigated carefully. Those points are 
addressed in the following paragraphs. 

4. 1 . Source morphology 

While we interpret and fit the observed flare as a single source, 
it cannot fully be excluded that the source is in reality a double- 
source in the 30-45 keV energy band as hinted by Pixon (see 
Sect. l3.1.2l and Fig.|4|, representing the footpoints of a flare loop. 
Since CLEAN finds indications of this structure, as well, it is un- 
likely that it is purely over-resolution caused by Pixon. However, 
visibility forward fitting was unsuccessful in fitting two sources 
reliably. Therefore, the single Gaussian model represents a good 
fit to the data but the spatial morphology in the third energy band 
is expected to have an influence on the fitted source size. As can 
be seen in Fig. [5] the size in the third energy band is somewhat 
larger than what would be expected from the general trend of in- 
creasing and decreasing size. While this might enhance the effect 
of increasing and decreasing source size, the underlying cause is 
likely to be albedo. At the same time, the size in the fourth en- 
ergy band (45-75 keV) is smaller than expected from any of the 
models. This could be caused by the smaller signal to noise ra- 
tio (weaker source) in this energy band that affects the result of 
visibility forward fitting (see Fig. [3}. This is also reflected in the 
uncertainty of the size in this energy band. 

In principle, the observed source could also be a very dense 
coronal source, as w as found in other RHESSI observations (e.g. 
IVeronig et al.ll2005ah . This scenario can, however, be ruled out 
as follows. The spectrum indicates a rather weak thermal com- 
ponent with emission measure of 7.4 x 10 46 cirT 3 . Assuming a 
true source FWHM of 7 arcsec (5.1 x 10 8 cm), the density of the 
source would be 2.8 x 10 I0 cirT 3 . This corresponds to a stopping 
depth of 3.2 x 10 9 cm for an electron with energy 30 keV which 
is much larger than the extent of the source. The electrons will 
therefore be stopped much lower down in the atmosphere, indi- 
cating that the HXR emission is indeed chromospheric and not 
coronal. 



4.2. Directivity and primary source size 

The FWHM of the observed source is found to increase and de- 
crease with energy with a peak between 30 ke V and 40 ke V. This 
is contr ary to what o ne would expect in the classical thick target 
model dBrownlll971l) . Collisional energy losses would lead to a 
smaller extent of the source along the direction of propagation 
(parallel to the magnetic field) as a function of energy. Since 
electrons with higher energies penetrate deeper into the target, 
the extent of the source perpendicular to the magnetic field is ex- 
pected to decrease as a function of energy in a converging mag- 
netic loop. An overall decrease of the source size with energy 
would therefore be exp ected. Such behavio ur was observed by 
iKontar et all (120081) and lKontar et al l (120101) in a limb event for 
which albedo has little effect. 

The most obvious explanation for the observed behaviour 
presented in this paper is albedo. The contribution of backscat- 
tered flux is negligible below 10 keV due to photoelectric ab- 
sorption. Photons with energies larger than 100 keV penetrate 
deeply into the atmosphere and are lost to the observer. The 
albedo contribution is therefore largest between 10 keV and 100 
keV (the energy range observed here) with a maximum around 
30 - 40 keV - where the maximum in source size is measured. 
The observed increase and decrease of the size is therefore con- 
sistent with the expected increase and decrease of the backscat- 



tered component. The observed size depends on two main pa- 
rameters, the true source size, as well as the directivity. A third 
quantity that has to be considered is the height of the source 
above the photosphere. In the simulations we assumed a height 
of 1 Mm. This is a reasonable assumption in a standard at- 
mosphere. Values of this order were also found observationally 
(Saint-Hilaire et ai] l2010l: Hpntar et al l2010l: lAschwanden et al.l 
2002). Since we cannot directly fit the true source size nor the 
directivity, we infer one assuming the other and compare the 
resulting size behaviour to the observations. As shown in Fig. 
[5] the data are consistent with a true source size of 7 arcsec 
for a downward directivity o f a = 2, a value that is consistent 
with previous observa t ions (iKand Il974t IVestrand et all 11987b 
iKasparova et al.l 120071: IKontar & Brownl 120061) . On the other 



hand, the inferred true source size is rather lar ge for a footpoint, 
compared with previous measurements (e.g. Dennis & Pernak 
2009). Even for a large directivity such as a = 5 the true source 
size could be of the order of 6 arcsec. It is conceivable that the 
true source is smaller if the velocity distribution of the X-ray pro- 
ducing electrons is not beam-shaped but rather takes the form of 
a "pancake" distribution. In this case, the downward directivity 
would be the same, but with a higher contribution of reflected 
photon flux from the sides, leading to a larger albedo patch and 
therefore larger source size. 



4.3. Direct and indirect imaging 

We chose an indirect approach to measure the albedo because 
the albedo patch is faint compared to the main source and is 
difficult to image directly. In CLEAN, any albedo component 
will most likely be masked by side-lobes of the main source. It 
has been suggested that it might be possible to image the albedo 
source directly using Pixon. We investigated this possibility us- 
ing simulated flare data. Starting from simulated photon maps 
of a Gaussian source with albedo (no background), calibrated 
event lists were computed. The RHESSI software was then used 
to reconstruct Pixon images from the calibrated event lists. For a 
small primary source (2 arcsec) and very high countrates, Pixon 
seems to image the albedo component in the form of "wings" 
at larger distances from the source centre. However, for weaker 
and larger sources (comparable with the data presented here) this 
does not seem to work. A primary source with 8 arcsec FWHM 
was simulated both with and without albedo component. Figure 
[6] illustrates the radial profile of the simulated source (summed 
over the y-axis) compared to the radial profile of the recon- 
structed Pixon image for the case of isotropic albedo (a = 1, 
left) and without albedo (i.e. pure Gaussian a — 0, right). The 
Pixon image depends strongly on certain parameters used for 
the image reconstruction, in particular the Pixon sensitivity. Two 
main points should be noted. First, in the albedo case Pixon finds 
"wings" at larger distances from the source centre that can be 
interpreted as the albedo component. However, in the a = sce- 
nario (no albedo), Pixon indicates some emission at larger dis- 
tances from the source centre that was not existent in the initial 
simulated image and must be artefacts from the image reproduc- 
tion. In real observations it will therefore be difficult to judge 
whether such "wings" are truly albedo or just imaging artefacts. 
Secondly, as has been known previously, Pixon has a tendency 
to over-resolve sources. This is clearly demonstrated here where 
a non-existent substructure is found in a source that is a single 
Gaussian. Another possible approach is to fit the albedo patch 
directly using visibility forward fitting, i.e. to fit two sources one 
of which is the (small) main source, the other the (large) albedo 
patch. In one of the simulations it was possible to fit one smaller 
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Fig. 6. Radial profiles (summed over the y-direction) of simulated photon maps including albedo (directivity a = 1, left) and without 
albedo (right). The thick blue line indicates the radial profile through the Pixon images computed from the simulations. The inset 
shows contour plots of the Pixon images. 



source (the size of the initial Gaussian) plus a larger Gaussian 
source (~ 15 arcsec FWHM) as the albedo component. Again, 
this is dependent on good count statistics and was not success- 
fully applied to the observed flare data. 



5. Conclusions 

In this work we discuss a new method to measure albedo in 
imaging. Since direct imaging of the faint albedo source is dif- 
ficult, an indirect approach via imaging spectroscopy of the 
size (second moment of the X-ray source distribution) is bet- 
ter suited. The method combines the energy dependence of the 
albedo with the ability of visibility forward fitting to measure 
large scale sources structures. This is a step forward compared 
to previous studies which only used o ne energy band to image 
albedo (eg. ISchmahl & Hurfo rd 2002). We emphasise that dif- 
ferent energy bands have to be analysed simultaneously to dis- 
tinguish potential extended primary sources from albedo. 

The observed size versus energy dependence contradicts the 
expected behaviour in a classical thick target model. It can most 
easily be attributed to albedo. Given the albedo interpretation of 
the source, we find that the true source size is likely to be of the 
order of 7 arcsec for realistic values of the directivity in the order 
of a — 2 which would be co nsistent with earlier observational 
studies of anistrppy of albedo (iKand i974tlVestrand et all l987: 
Kasn arova et al.1 l2007t iKontar et al.1 120061) . A model involving 
increasing directivity with energy can be ruled out by the data. 

Although the proposed method is limited to simple source 
morphologies at the moment, it clearly demonstrates the influ- 
ence albedo has on observed X-ray sources and the derived sizes, 
a fact that has to be considered when inferring physical proper- 
ties of flares. Further, it opens the way to determining the di- 
rectivity of X-rays and inferring the directionality of accelerated 
electrons not only from spectra, but from imaging, for the first 
time. If the uncertainties could be reduced by a factor of 3, a con- 
fident rejection of a wider class of beam transport models would 
be possible. This might be achieved in a (future) flare with higher 



count rates or/and by further improving RHESSI visibility anal- 
ysis, as well as in measurements made by future instruments. 
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